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ABSTRACT 

A method utilizing vector representation is investigated for deter- 
mining a faulty element in passive and active networks by simple 
external measurements. 

A large system may be considered as an interconnection of a number 
of subnetworks. By utilizing the relationships between the magnitudes 
of a transfer function at various frequencies and the deviations of a 
circuit element, the fault simulation curves can be drawn. The fault 
identification regions are defined from the fault simulation curves. 

A fault identification matrix is constructed corresponding to the de- 
fined fault identification regions. 

The fault identification matrix, when premultiplied by a vector 
whose components are measured from a netivork, yields another vector 
whose components identify a network element which is faulty. 

A test procedure for the fault identification method is presented 
and verified by experiments. 
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I. INTRODUCTION 



Recent advances in the design of electronic systems have resulted in 
an ever-increasing complexity in such systems. Despite the advances in 
component reliability and preventive maintenance, it appeared that a 
great many unpredictable or random failures do occur. Such complexity 
in modern equipments used by the industry and the military has resulted 
in a serious problem of providing adequate maintenance. One of the most 
difficult tasks is the location of malfunctions v;hich, until recently, 
has been left entirely in the hands of technicians. 

Some teichniques for fault identification have been studied by a 
number of investigators [1,2,3,5,6,7,8,9,10] and a computer program for 
fault simulation of passive linear networks without mutual inductance 
has been given -in Ref. [2]. 

In this paper, the problem of fault identification, which require-s 
the least number of external test points of the system, is investigated 
and a procedure which will automatically identify a fault or faults in 
passive and active linear networks is presented. Toward achieving this 
ultimate goal, the sensitivity of a network function due to element devi- 
ations is ‘studied. The method, which requires matrix operations, is 
introduced. In this method, the fault identification regions at a set 
of selected test frequencies are defined depending on the network sensi- 
tivity and the fault identification matrix is generated. A vector is 
then obtained by premultiplying the fault identification matrix by 
another vector, and the resulting vector identifies the fault. 

The single element fault identification method is discussed and the 
combined element fault identification method is introduced by utilizing 
the single element fault identification method. 
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The computer programs, which identify a fault element in a properly 
partitioned network, using the proposed method are written in FORTRAN IV 
and included in Appendix A. The implementation of the computer programs 
is annexed in Appendix B. 
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II. PRELIMINARY CONSIDERATIONS 



A. DEFINHION OF A VECTOR 

Consider a two-port network with the' external terminals and variables 
as defined in Fig. 2-1. The forward voltage ratio transfer function, 



1 ' o 

Fig. 2-1. 

2-port netv/ork. 

T12, is defined as the ratio of the output voltage to the input voltage 
with the output current set to zero. 




T12 




(^-1) 



The magnitude of T12, designated by ltl2| , is given by either a 
function of the frequency with fixed element values or a function of an 
element value with a fixed frequency. That is, at a set of n frequen- 
cies, with fixed element values, a set of n different magnitudes of a 
transfer function, T12, can be obtained. Similarly, corresponding to a 
set of k deviations of an element value, with a fixed frequency, a set 
of k different magnitudes of T12 can be determined. 

Let the network contain only the following types of passive elements: 
resistors, inductors and capacitors. Consider an element v/hich has a 
nominal element value E with a deviation designated by aE. , where the 

J 

subscript j implies the j-th deviation. 
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With a change in element value aE. of E, n values of |T12| can 

J 

be obtained corresponding to n different frequencies. These values of 
|T12| are designated by 

|T12|., , IT12I.2 . - , |T12|j, . - , |T12|j„ 

where the first subscript indicates the j-th deviation and the second 
subscript indicates the i-th frequency. Next, with a nominal element 
value of E, n values of |T12| can be obtained corresponding to n dif- 
ferent frequencies and designated by 

|T12|j, , |T12|(,2 . — , |T12|(,, . — , |T12|(,„ 

The normalized value of |T12| . . is defined as 

' 'JT 

|T12|j, 

where i = 1, 2, — , n. 

For example, with a deviation aE. of E and a set of 5 frequencies, the 

J 

5 values of the normalized |T12|,. are 



With another 
the 5 values 



|T12|„ |T12|i2 |T12|,5 

■ iTWi- ■ ’ 

deviation, aE^ , of E and the same set of 5 frequencies, 
of the normalized |T 12 | 2 .j are 



IT12I2, 1T12I22 ■ IT12I25 

rnri^ ■ TniT^ ’ ’ 
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Let a vector X be denoted by 



X = ( x^, x^, — - , x^. , , x^, c )' (2-2) 

where X is an array of the n normalized values of |T12) with respect 
to a set of n frequencies (with a deviation of aE.) and a constant c: 
x^- denotes j-j-i g (i = 1,2, — , n) and c is unity; 

X is n dimensional because the last component is always unity; 
the superscript (') means a transpose of an array of T or a matrix. 

B. CONSTRUCTION OF FAULT IDENTIFICATION MATRIX AND 
GENERATION OF VECTOR-Z 

Consider a simple two-port network that contains a resistor, R, and 
a capacitor, C, as shown in Fig. 2-2. Let the nominal value of R be 
1 ohm and that of C be 1 farad. The value of C can be changed from 



r\ 




r 


n _ _ _ 


o 

II 




( 






V. 



Fig. 2-2. 

2-port RC network. 



zero to infinity and 3 nonzero finite values of C are indicated as shown 
in Fig. 2-3. 

0.0 0.5 1.2 1.5 

I 1 1 1 5 



Fig. 2-3. 
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In general, consider that the value of a parameter can vary from 
zero to infinity and this range of values can be divided into a number 
of subdivisions, each of v/hich shall be called the j-th deviation of the 
parameter value, v-/here j = 1, 2, — . For example, the value of C 
(betv/een 0 to “) is divided into four subdivisons, that is, the 1-st 
deviation is 0 £ C < 0.5 , the 2-nd deviation is 0.5 < 1.2 and 

similarly, 1.2 < 1.5 and 1.5 represent, respectively, the 

3-rd and 4-th deviation. 

The normalized |T12| is given by 



The values of the normalized |T12| corresponding to the 1-st deviation 
through the 4-th deviation are shown in Fig. 2-4 according to a set of 
3 test frequencies. The values of the normalized |T12| corresponding 
to the j-th deviation where j = 1,2, 3, 4 are best seen by a coordinate 
plane whose ordinate represents a frequency axis and abscissa represents 
a normalized |T12| axis as shown in Fig. 2-5, 

Let a "region-j", designated by R-j , represent an area on the fre- 
quency-normalized |T12| coordinate plane which corresponds to the j-th 
deviation of an element. For example, R-1 represents the area to the 
right of line A-0 (Fig. 2-5) on the frequency-normalized |T12| coordi- 
nate plane and corresponds to the 1-st deviation of C. R-2 represents 
an area between line B-0 and line A-0, including line A-0, on the fre- 
quency-normalized |T12| coordinate plane and corresponds to the 2-nd 
deviation of C as shown in the figure. 
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1 -St 

devi ation 



2-nd 

deviation 



C f- 



w=l 



w=2 



w=4 



0. 


5 

1 




1.265 




1.582 




1.845 







3-rd 

deviation 



4-th 

deviation 



1.2 



1 .5 



1 


0.906 j 


1 


T 


0.861 




2 






3 


> 


0.845 > 



2 ^ 



'3 ^ 



0.785 



'1 ^ 



0.707 



0.68 



3 > 



Fig. 2-4. 

The normalized |T12| corresponding to the j-th deviation. 




Fig. 2-5. 



The frequency-normalized |T12j coordinate plane. 
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A "normal region" is defined as a region on the frequency-normal ized 
|T12| coordinate plane where a corresponding element deviation is in a 
predefined normal range of deviation. The region outside of this normal 
region is defined as a "faulty region". For example, if the 2-nd devi- 
ation is a predefined normal range of deviation, R-2 is a normal region. 
The regions outside of R-2, that is, R-1 , R-3 and R-4 are considered to 
be faulty regions. 

Let be a variable in a region R-j of domain x. If the value of 
x^ is in a normal region, the corresponding element is considered to be 
normal. For example, let R-2 be a normal region. If the values of x^ 
where i = 1,2,3 is in R-2, the element C is’ considered to be normal. 
Otherwise the element with the values represented by x^. is considered 
to be faulty. 

Let a "boundary point value", abbreviated by B.P.V., represent a 
value of the normalized |T12| which lies between two regions. The 
boundary point values of a region may or may not belong to that region. 
For example, the boundary point values between R-1 and R-2 lie on line 
A-0 and these values do not belong to R-1 but belong to R-2 as shown 
in Fig. 2-'5. The 3 boundary point values corresponding to x.j , x^ and 
x^ in R-1 are 1.265 , 1.582 and 1.845 and belong to R-2. 

A "fault identification matrix", designated by M, is defined as a 
matrix whose number of columns represent the number of regions and the 
number of rows represent the number of components of X as shov/n in 
Fig. 2-6. 

A vector " Z" " is defined as 



V = X'M 



(2-4) 
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M = 



(No. of regions 
R-1 R-2 — 



No. of columns) 

R-j — - R-r 
• _x_ I 



(No. of X components = No. of rov/s) 
Fig. 2-6. 

Fault identification matrix. 



1 

2 

i 

n 

n+1 



A method for the construction of M and the generation of T in the 
n-dimensional domain x with r-regions is presented as follows. The 

r-regions R-1, R-2, , R-j, , R-r and p boundary point values 

b.j , , bp , assigned by positive numbers, are defined in Fig. 

2-7. 

Corresponding to the boundary point values, the r-regions can be 
arranged with respect to the vector component x^. as shown in Table 2-1. 
Let the regions, arranged as shown in Table 2-1, be called the "identi- 
fication regions". 

A fault identification matrix, M, will be constructed corresponding 
to the defined identification regions. For every B.P.V. associated 
with each vector component, x^. , a sub-fault identification matrix is 
formed. These sub-fault identification matrices, designated by m-j , m 2 
, have the same dimension as H. 
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R-1 



R-2 



I 

I 

I 



R-j 



I 

I 

I 



R-r 




Fig. 2-7, • 

n-dimensional domain x with r-regions 




Identification regions. 
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The following rules are given for the construction of a fault identi- 
fication matrix. 

Rule 1 : 

To construct a sub-fault identification matrix for a particular 
x^ , the following procedures are used. 

(a) . If the value of x^. in R-j is greater than the selected 
B.P.V., that is associated with x^ , put +1 in the i-th row of j-th column 
of the sub- fault identification matrix. 

(b) . If the value of x^. in R-j is less, than the selected 
B.P.V. , that is associated with x^. , put -1 in the i-th row of j-th column 
of the sub- fault identification matrix. 

(c) . To generate the entry of the j-th column of the last 
row, multiply the selected B.P.V. by the entry value of the j-th column 
of the i-th row and change its sign. 

(d) . All entries other than those in the i-th row and the 
last row are set to zero. 

Rule 2: 

The fault identification matrix, M, is constructed as a matrix 
sum of al 1 • sub-faul t identification matrices. 

Using rules 1 and 2, the sub-fault identification matrices and 
M for the identification regions in Table 2-1 are constructed as follows. 
With respect to x.j , the sub-fault identification matrices are: 
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for B.P.V. of b 



(p-1)’ 



m 



1 



R-1 R-2 

- 1 . - 1 . 

0 . 0 . 



0 . 0 . 



R-j - — 



R-r 



+ 1 . 

0 . 



1 

2 



0 . 



i 

n 



■'’(p-i) 



n+1 



and for B.P.V. of 






R-1 


R-2 


+1. 


-1. 


0. 


0. 



R-j 



0 . 0 . 




+b 



2 



R-r 

+ 1 . 



0 . 



0 . 




1 

2 

i 

n 

n+1 



With respect to x^, the sub-fault identification matrices are: 
for B.P.V.- of b-, 



R-1 


R-2 


— R-j — 


R-r 


0. 


0. 




0. 


-1. 


-1. 




+1. 


0. 


0. 




0. 



1 

2 

i 

n 



+b. 



+b. 



n+1 
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and for B.P. V. of , 


R-1 


R-2 


— R-j — 


R-r 






0. 


0. 




0. 


1 




+1. 


+1. 




+1 . 


2 












i 




0. 


0. 




0. 


n 




.■'i 


-t-i 




■”1. 


n+1 


Finally, with respect to 


x^, the 


sub-faul t 


identification 


matrix for 


B.P.V. of bp is 


R-1 


R-2 


— R-j — 


R-r 




= 


0. 


0. 




0. 


1 




0. 


0. 




0. 


2 












i 




-1. 


-1. 




+1. 


n 




/'’p 






■■’p. 


n+1 


where means the last : 


sub-fault 


identification matrix. 






The resulting M is given by 










M = m^ 


+ m 2 + m 


3 + m4 + - 


+ m^ 








R-1 


R-2 


— R-j — 


R-r 




M = 


0. 


-2. 




+2. ~ 


1 




0. 


0. 




+2. 


2 












i 




-1. 


-1. 




+1. 


n 






A2 






n+1 
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■I II II II II 



where A, " bjp.,, - + 2b2 - b, e — + 

and \ ‘ - t’(p.i) - 2bg - b, + b^ . 

Next, the application of M for the generation of 1 is discussed. 

Let the normalized values of lT12| of a system be x-j , X 2 > , x^ . 

X is in the form of 

X “ E x-| , x£ j , > ^p’ J 

T' is generated as follows. 

T = X'M = [ x^, X2, , x^, 1. ] M 

- E ’ ~2x-|+ ”^n"*^2 ’ ’ 2x-|+2x2 ~~ "*”^n "^r ^ 

~ E Z-| , Z 2 , , Zj , - , z^ 3 

Z has the same number of components as the number of regions and Zj 
serves as an indicator of the region in v/hich the system is operating as 
given by the following rule. 

Rule 3: 

If a system is operating in R- j , the component z. of 2 is the alge- 

. • vj 

braically largest component in Z. 

The construction of M and the generation of Z are illustrated in 
the following examples. 

Example 2.1 . Construct M and generate Z in the following one 
dimensional domain x with three regions. 



R-2 



X < 5.0 



R-1 



R-3 



5.0 



5.0 < X < 10.0 



10.0 < X 
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Corresponding to the given boundary point values, the identification 
regions are as follows. 

R-1 R-2 R-3 

X 5.0j<x<10.0 x<5.0 lO.O^x 

Using rules 1 and 2, the sub-fault identification matrices and M for the 
identification regions are constructed. 

The sub-fault identification matrix for B.P.V. of 5.0 is 





R-1 


R-2 


R-3 


m^ = 


“ +1. 


-1. 


• +1. ' 




-5.0 


+5.0 


-5.0 



and the sub-fault identification matrix for B.P.V. of 10.0 is 









R-1 


R-2 


R-3 




m2 


= 


"-1. 


-1. 


+1. ' 








_+10.0 


+10.0 


-10. 0_ 


resulting 


M is 


given 


by 












R-1 


R-2 


R-3 


M = 


+ m2 


= 


0. 


-2. 


2. " 








5.0 


15.0 


-15. 0_ 



Suppose that the measured value of x is 13.0. X is in the form of 

X ' = E X , 1 . ] 

[ 13.0 1. ] 
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Thus, V is obtained 



Z' = X'M = [ 


13.0 


1. ] 


~0. 


-2. 2. ' 








_5.0 


15.0 -15. 0_ 


= [ 


5. 


-11. 


11. ] 




[ 


^.1 ’ 




^3 ^ 





in which is the algebraically largest component. According to rule 3, 
the system is operating in R-3, The measured value, of x, indeed, belongs 
to R-3. 



If the measured value of x is 6.0, then 



V = 



= X'M = [ 6.0 


1. ] 


1. 


-2. 


2. 






. 5.0 


15.0 


-15.0_ 


= [ 5. 


3. 


1 1 

CO 





= r 2l . 



•2 ’ 



^3 ^ 



The fact that z-j = 5 , the largest of all the entries in Z' 
that the system is operating in R-1. . ' 



indicates 
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Example 2.2. Construct M and generate Z in the following three 



dimensional domain x with four regions. 



R-1 



R-2 



R-3 



R-4 



1.5 2.0 2.5 



K- X, 



^2 ^ 



4.0 5.0 



Xi > 



X, 



« — X, 



X. 



The boundary point values are 1.5 , 2.0 , 2.5 , 4.0 and 5.0 . The 
identification regions are defined as follows. 





R-1 


R-2 


R-3 


R-4 


^1 


2.0 < x-| < 4.0 


x.| < 2.0 


4.0 < x.| 


2.0 < x.| < 4.0 


^2 


1.5 < X 2 < 2.5 


1.5 < X 2 < 2.5 


2.5 < X 2 


1.5 < X 2 < 2.5 




x^ < 5.0 


x^ < 5.0 


X 3 < 6.0 


5.0 < x^ 
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V/ith respect to , the sub-fault identification matrices are: 
for B.P.V. of 4.0, 



R-1 


R-2 


R-3 


R-4 




-1. 


-1. 


1. 


-1. ' 


1 


0. 


0. 


0. 


0. 


2 


0. 


0. 


0. 


0. 


3 


4. 


4. 


-4. 


4. 


4 



and for B.P.V. of 2.0, 



R-1 


R-2 


. R-3 


R-4 




1. 


-1. 


1. 


1. " 


1 


0. 


0. 


0. 


0. 


2 


0. 


0. 


0. 


0. 


3 


-2. 


2. 


-2. 


-2. 


4 



With respect to X 2 > the sub-fault identification matrices are: 
for B.P.V. of 2.5, 



R-1 


R-2 


R-3 


• R-4 




0. 


0. 


0. 


0. 


1 


-1. 


-1. 


1. 


-1. 


2 


0. 


0. 


0. 


0. 


3 


2.5 


2.5 


-2.5 


2.5_ 


4 



22 



and for B.P.V. of 1 .5 





R-1 


R-2 


R-3 


R-4 






% = 


0. 


0. 


0. 


0. “ 


1 






1. 


1. 


1 . 


1. 


2 






0. 


0. 


0. 


0. 


3 








-1.5 


-1.5 


-1.5^ 


4 




With respect to x^, the sub-fault identification matrix 


for B 


.P.V. of 




5.0 is 
















R-1 


R-2 


R-3 


R-4 






"^5 = 


0. 


0. 


0. 


0. 


1 






0. 


0. 


0. 


0. 


2 






-1. 


-1. 


-1. 


1. 


3 






5. 


5. 


5. 


-5. _ 


4 




The resulting M is given by 
















R-1 


• R-2 


R-3 


R-4 




M = m^ + m2 + m^ + m^ 


1 + = 


■ 0. 


-2. 


2. 


0. ~ 


1 






0. 


0- 


2. 


0. 


2 






-1. 


-1. 


-1 . 


1. 


3 






_ 8. 


12. 


-5. 


-2._ 


4 


If the measured components of X are 


^1 " 


3.0 , X 2 = 2. 


.0 and X 2 = 3.0 



X' becomes 



V 



X'M 



[ 3. 



E 3. 2. 3. 


1 . ] 


1 1 

II 

1 1 

cn 


3 . 2 . 1 . ] 


N 

1 — 1 
11 


> ^2 > 3 
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The component z-j is the albebraical ly largest component. According to 
rule 3, the system is operating in R-1 . Similarly, if the measured 
components of X are x.j = 3.0 , ^2 = 2.0 and = 6.0 , Z is given by 

r = X'M = [ 3. 2. 6. 1 . ] M = [ 2. 0. -1. 4. ] 

I 2 i , , 23 , 24 ] 



The component z^ indicates that the system is operating in R-4. 
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III. THE PROPOSED METHOD 



A. PURPOSE 

The purpose of the proposed method is to investigate a feasible method 
to determine whether fault detection and isolation can be accomplished by 
measuring only the magnitude of the system transfer frunction. In this 
method, the least number of test frequencies and test points, such as the 
normal input and output terminals of the system, will be used. 

B. BASIS FOR THE PROPOSED METHOD 

Consider again the 2-port network shown in Fig. 2-1. The open circuit 
impedance functions are represented by the equations of the voltages as 
linear functions of the currents in matrix form as follows. 



1 




I 

N 


^12 




'h' 


1 

CM 




^21 


^22 




'2 



The network of the type being considered can be externally character- 
ized by a transfer function in terms of open circuit impedances. It is 
easily found that the forv;ard voltage ratio transfer function, T12, and 
the reverse voltage ratio transfer function, T21, (defined as the ratio 
of the input voltage to the output voltage with the input current set to 
zero) are 

T12 = ^ (3-2) 

^11 

T21 = ^ (3-3) 

^22 
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In general, a transfer function can be expressed as a rational func- 
tion of the complex frequency, s, and a network element, E, as follows. 



ra 

K n (s + z^. ) 

T(s,E) = (3-4) 

n 

n (s + p.) 
j=i J 



where zero, z. , and pole, p. , and K are functions of E and m and n are 
the number of zeroes and poles respectively. 

Taking the logarithm of Eq. (3-4), the following equation is obtained 

m n 

In T(s,E) = In K + Z ln(s + zj - Z ln(s + p.) (3-5) 

i=l ^ j=l J 

And taking the derivative of Eq. (3-5), the following equation is 
obtained. 



d(ln T(s,E)) 



d(T(s,E)) 

T(s,E) 



= i ^ 

K dE 



m 1 dz. 

Z — 

i = l s+z^. dE 



n 1 dp. 

Z — ' — — dE (3-6) 
j=l' s+p. dE 

J 



A "transfer function sensitivity" is represented by Eq. (3-6). The 

transfer function sensitivity is closely related to pole and zero sensi- 
dp. dz. 

tivities where 3 ^d defined as pole sensitivity and 

zero sensitivity respectively. In other words, the variation of T(s,E) 
may take the form of a change in location of a pole and/or a zero of 
T(s,E). 
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From Eq. (3-4) and Eq. (3-6), it can be seen that a change in any 
one of the pole-zero configurations will cause a change in one or possi- 
bly both of the voltage ratio transfer functions. Since T12(s,E) and 
T21(s,E) depend on the value of E, it is not hard to consider how a 
change of an internal element is related to a change of an external 
voltage ratio transfer function. 

This is the basic philosophy of the proposed method for fault 
identification. 



C. PROPOSED TEST PROCEDURE 

A fundamental fault testing procedure for the proposed method is 
given as follows. 

Step (1). Partition the system. 

Step (2). Generate a symbolic system function. 

Step (3). Select a set of test frequencies. 

Step (4). Simulate a fault. 

Step (5). Determine the fault identification regions. 

Step (6). Construct the fault identification matrix. 

Step (7). Generate X. 

Step (8). Generator. 

Step (9). Identify the fault. 

(a) . Single element. 

(b) . Combined element. 

Step (2), step (3) and part of step (4) have been studied and a 
computer program of these steps for the passive linear network without 
mutual inductance was presented by J. D. Courville in Ref. [2]. The 
detailed discussion of these steps will be omitted. The other steps are 
discussed in detail in the next chapter illustrative examples. 
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IV. NETWORK FAULT IDENTIFICATION 



A. NETWORK CONFIGURATIONS AND TRANSFER FUNCTIONS 

As illustrated in the foregoing discussions, there is a certain 
relationship between the transfer function and a network element devia- 
tion. 

To illustrate the proposed method, two examples for passive network 
configurations are presented. 

Exampl e 4. 1 : 




2 

2 ' 



1/Eal = R.| = 2 ohms 
1/Ea2 = L =2 henry 
Ea3 = C '= 3 farads 
1/Ea4 = R^ = 4 ohms 



Fig. 4-l(a). 
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Example 4,2: 




Eb3 = C =1 farad 
1/Eb4 = L 2 = 1 henry 



Fig. 4-1(b). 



For example 4.1, the transfer function is 



T12- 



R-jLCs + (R-|R 2 C+L)s + R.| + R 2 



1 

3s^ + 6.5s + 1.5 



(s+0.262)(s+1.90) 
where K.| is 1/3 . 



(4-la) 
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For example 4.2, the transfer function is 



T12 = 



L^s 

L^L^Cs^ + RL^Cs^ + (L.|+L2 )s + R 



3 2 

s'^ + s^ + 2s + 1 



K2S 



(s+0.57)(s^+0.43s+1.7623) 



(4-lb) 



where K 2 is 1 . 

The network semi logarithmic plots for the magnitude characteristic 
of the transfer function are shown in Fig. 4-2a for Eq. 4-la and Fig. 

4-2b for Eq. 4-lb. 

B. NECESSARY DATA 

All necessary data for the proposed method are based on a set of 
chosen test frequencies. 

A "fault simulation curve" is defined as a response' curve of the 
normalized |T12| variations due to the simulated deviations of an 
element with respect to a chosen test frequency. The fault simulated 
data for lT12| can be obtained from the computer program in Ref. [2]. 

The fault simulated data for lT12| and the normalized lT12| are 
tabulated in Tables 4-la and 4-lb for Eq. 4-la and Eq. 4-lb respectively. 

The fault simulation curves are drawn in Graph 4-1 a-1 thru Graph 
4-la-3 for example 4.1 and Graph 4-lb-l thru Graph 4-lb-3 for example 4.2. 
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.symptotic magnitude versus frequency plot for 
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c. DETERMINATION OF FAULT IDENTIFICATION REGIONS 



The proper determination of the identification regions is the most 
laborious part of the proposed method. 

For example 4.1, a set of 3 test frequencies is chosen as follov/s. 

= 2.00 radians/sec 

W 2 = 1.02 radians/sec 

= 0.02 radians/sec 

For example 4.2, a set of 3 test frequencies is chosen as follov/s 

= 3.20 radians/sec 

W 2 = 0.70 radians/sec 

= 0.05 radians/sec 

For a test frequency, there are eight basic identification regions, 
four normal regions and four faulty regions, because four elements are 
contained in the examples. That is, each element has two regions, a 
normal region and a faulty region. The four normal regions are not the 
concern of the propsed method. The four faulty regions may be sub- 
divided into many faulty regions to identify a fault. These subdivisions 
of the faulty regions depend on the transfer function sensitivity. The 
subdivided faulty regions are called the fault identification regions. 

The normal regions for the networks are defined as in Table 4-2a for 
example 4.1 and Table 4-2b for example 4.2 according to the fault simu- 
lation curves. The regions outside of the normal regions are considered 
to be faulty regions. The fault identification regions are defined as 
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DEV 


MAG-Tl 


MAG-T2 


MAG-T3 


NOR-Tl 


N0R-T2 


N0R-T3 


(percent) 


(w=2.0) 


(w=1.02) 


(w=0.02) 








1 


90 


0. 11247 


0.26665 


0.79074 


1.8795 


1.8199 


1.1896 




70 


0. 10C92 


0.24123 


0.77166 


1.6865 


1 .6464 


1 .1609 




50 


0. C8928 


0.21504 


0.74875 


1 . 4920 


1 .4677 


1.1264 




30 


0. C7756 


0.18812 


0.72074 


1.2961 


1 .2840 


1. 0843 




10 


0. C6576 


0. 16054 


0. 68572 


1.0990 


1 .0957 


1.0316 




-10 


0. G5390 


0.13235 


0.64069 


0. 9008 


0.9033 


0. 9639 




-30 


0.04199 


0.10262 


0.58067 


0. 7017 


0.7072 


0.8736 




-50 


0.03003 


0.07442 


0.49672 


0.5018 


0.5079 


0.7473 




-70 


0. C1603 


0.04485 


0.37120 


0. 3014 


0.3061 


0.5584 




-90 


C. 00602 


0.01500 


0. 16370 


0. 1005 


0.1024 


0.2463 


2 


90 


0.07451 


0.15650 


0.66464 


1.2452 


1.0681 


0.9999 




70 


0.07267 


0.15555 


0. 66465 


1.2144 


1.0617 


0.9999 




50 


C.C7C24 


0.15418 


0. 66466 


1. 1737 


1.0523 


0.9999 




30 


0. C6699 


0.15211 


0.66467 


1.1195 


1 .0382 


0.9999 




10 


0. C6260 


0.14885 


0.66470 


1.0461 


1 .0160 


1.0000 




-10 


0.05662 


0.14247 


0.66473 


0. 9462 


0.9792 


1.0000 




-30 


0. 04851 


0.12404 


0.66477 


0.8106 


0.9149 


1.0001 




-50 


0. 03774 


0.11667 


0 . 66484 


0. 6306 


0.7963 


1.0002 




-70 


0. 02413 


0.08429 


0.66497 


0.4032 


0.5760 


1.0004 




-90 


0. CC830 


0.02149 


0.66482 


. 0.1386 


0.2149 


1.0002 


3 


90 


0.03131 


0.07739 


0.65942 


0. 5232 


0.5282 


0.9920 




70 


0. C5502 


0.08652 


0.66088 


0. 5853 


0.5905 


0.9942 




50 


0.C3974 


0. 09805 


0.66217 


0.6641 


0.6692 


0.9962 




30 


0.04591 


0.11208 


0.66331 


0. 7673 


0.7718 


0.9979 




10 


0. C5435 


C. 12241 


0.66426 


0. 9082 


0.9106 


0.9994 




-10 


0. C6656 


0.16238 


0.66509 


1.1122 


1 .1082 


1.0006 




-30 


0.08576 


0. 20630 


0.66574 


1.4331 


1.4094 


1.0015 




-50 


0.12017 


0.28008 


0.66621 


2.0081 


1.9116 


1.0023 




-70 


0. 19776 


0.41447 


0.66651 


3. 3 047 


2.8288 


1.0027 




-90 


C. 45038 


0.61200 


0.66665 


7.5262 


4.1770 


1.0029 


4 


90 


0. C3991 


0.12301 


0. 51212 


0. 6669 


0.8395 


0.7704 




70 


0. C4336 


0.12829 


0.53966 


0. 7246 


0.8756 


0.8119 




50 


0. C4733 


0.13364 


0.57032 


■ 0.7910 


0.9121 


0.8580 




30 


0.C5188 


0.13894 


0.60467 


0. 8670 


0.9483 


0.9097 




10 


0. C5704 


0.14407 


0.64342 


0.9532 


0.9833 


0.9680 




-10 


0. 06277 


0.14866 


0.68746 


1.0489 


1.0160 


1.0342 




-30 


0.06884 


0.15213 


0.73795 


1.1504 


1.0451 


1.1102 




-50 • 


0.07478 


0.15668 


0.79642 


1. 2497 


1-.0694 


1.1981 




-70 


0. 07973 


0.15932 


0.86491 


1.3324 


1.0874 


1.3012 




-90 


0.08264 


0.16091 


0.94624 


1.3809 


1 .0982 


1.4235 



Table 4-1a. 

The magnitudes and the normalized magnitudes of 

T12 due to each element deviations for Example 4.1. 

( The nominal magnitudes are 0.05984 at w=2.0 , 
0.1465 at w=1.02 and 0.6647 at w=0.02 ) 
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Ebl 



Eb2 



Eb3 



Eb4 



DEV 


MAG-Tl 


MAG-T2 


MAG-T3 


NCR-Tl 


N0R-T2 


N0R-T3 


(percent) 


(w=3.2) 


(w=0 . 7 ) 


(w=0.05) 








90 


0.24394 


0.60490 


0.04991 


2. 1299 


1 .0142 


1.0009 


70 


0.21254 


0.60589 


0.04990 


1.8558 


1 .0158 


1.0008 


50 


C. 18275 


0.6C623 


0. 04990 


1.5957 


1 .0 164 


1.0007 


30 


0. 15445 


0.60512 


0.04989 


1 . 3485 


1.0146 


1.0005 


10 


0. 12751 


0.6C085 


0.04967 


1.1134 


1.0074 


1.0002 


-10 


0. 10186 


0.58945 


0 . 04985 


0. 8893 


0.9883 


0.9997 


-30 


0. 07739 


0 .56150 


0.04980 


0. 6757 


0.9414 


0.9987 


-50 


0. 05403 


0.49499 


0.04968 


0.4717 


0.829Q 


0.9963 


-70 


0.C3170 


0.35117 


0.04925 


0. 2768 


0.5888 


0.9877 


-90 


0.01034 


0.121C4 


0. 04468 


0.0903 


0.2029 


0.8961 


90 


0.11935 


0.64187 


0.09353 


1. 0420 


1.0762 


1.8759 


70 


0. 11886 


0.63708 


0.08398 


1.0378 


1.0681 


1.6843 


50 


C. 11818 


0.62C43 


0.07433 


1.0318 


1.0570 


1.4908 


30 


0.11718 


0 .62085 


0. 06460 


1.0231 


1.0409 


1.2956 


10 


0.11563 


0 . 60646 


0.05479 


1 . 0096 


1 .0168 


1.0989 


-10 


0.11309 


0.56365 


0.04492 


0. 9874 


0.9786 


0.9008 


-30 


0. 10852 


0.54525 


0.03499 


0. 9475 


0.9142 


0.7018 


-50 


0. 09938 


0.47653 


0.02502 


0. 8677 


0.7990 


0.5019 


-70 


C. C7692 


0.34967 


0.01503 


0. 6891 


0.5863 


0.3014 


-90 


0. C3530 


0.13439 


0.00501 . 


0.2908 


0.2253 


0.1005 


90 


0. C5440 


0.93146 


0. 04997 


0.4749 


1 .5617 


1.002? 


70 


0.06158 


0.83950 


0. 04995 


0. 5377 


1 .4075 


1. 0017 


50 


C.07C95 


0.75730 


0.04992 


0.6195 


1 .2697 


1.0012 


30 


0.08369 


0.66570 


0.04990 


0. 7307 


1.1497 


1.0007 


10 


0. 1C20C 


0.62398 


0.04987 


0. 8906 


1 .0462 


1.0002 


-10 


0.13056 


C.57C89 


0.04985 


1.1400 


0.9572 


0.9998 


-30 


0.18130 


0.52510 


0.04982 


1.5830 


0.8804 


0. 9993 


-50 


0.29628 


0 .48544 


0. 04980 


2.5869 


0.8139 


0.9988 


-70 


0.7985C 


C. 45090 


0. 04877 


6. 9719 


0.7560 


0.9782 


-90 


1. C2460 


0.42067 


0.04975 


8.9460 


0.7053 


0.9978 


90 


0. 10913 


0. 35960 


0. 02630 


0.9528 


0.6032 


0.5274 


70 


0.11115 


0. 39472 


0.02938 


0, 9705 


0.6618 


0.5893 


50 


0.11275 


0. 43707 


C. 03329 ■ 


0. 9844 


0.7328 


0.6677 


30 


0.11386 


0.48950 


0.03B40 


0.9942 


C.6207 


0.7701 


10 


0.11445 


0.55602 


0.04535 


0. 9993 


0.9322 


0.9095 


-10 


0.11448 


0.64307 


0.05537 


0. 9995 


1.0782 


1.1105 


-30 


0.11395 


0.76155 


0.C71C6 


0.9949 


1.2768 


1.4252 


-50 • 


0. 11289 


0.93135 


0.09911 


0. 9856 


1.5615 


1.9877 


-70 


0.11134 


1 . 19190 


0.16324 


0.9721 


1.9984 


3.2740 


-90 


0. 10935 


1.62772 


0.43908 


0. 9548 


2.7291 


8.8061 



Table 4-1b. 

The magnitudes and the normalized magnitudes of 

T12 due to each element deviations for Example 4.2. 

( The nominal magnitudes are 0.11453 at w=3.2 , 
0.5964 at w=0.7 and 0.04986 at w=0.05 ) 
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Fault simulation curve for Example 4.1 at w=2.0. 
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Fault simulation curve for Example 4.1 at w=1.02. 
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Graph 4-1a-3, 

Fault simulation curve for Example 4.1 at w=0.02. 
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Fault simulation curve for Example 4.2 at w=3.2. 
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Fault simulation curve for Example 4.2 at w=0.7. 
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Fault simulation curve for Example 4.2 at w=0.05. 
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Normal regions for Example 4.2. 
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Fault identification regions for Example 4.1. 



FAULT IDENTIFICATION REGIONS 
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Fault identification regions for Example 4.2. 



1n Table 4-3a for example 4.1 and Table 4-3b for example 4.2. In Tables 
4-3a and 4-3b, an arrow "f" means "more than" and an arrow "T" means 
"less than". 

D. CONSTRUCTION OF FAULT IDENTIFICATION MATRIX 

The sub-fault identification matrices and the fault identification 
matrix, M, for the fault identification regions defined in Table 4-3a, 
are constructed step by step according to the rules given previously. 

The boundary point values for x^ where i - 1,2,3 are tabulated in Table 
4-4. 





B.P.V.'s 


^1 


1.2 


1.1 


0.9 


0.8 


X 2 


1 .2 


1.05 


0.95 


0.8 


^3 


1.06 


1 .02 


0.98 


0.92 



Table 4-4. 



The sub-fault identification matrices with respect to x.j are repre- 
sented by m.| , m 2 , m^ and m^ , the sub-fault identification matrices with 
respect to X 2 are represented by m^, mg, m^ and mg and the sub-fault 
identification matrices with respect to Xg are represented by mg, m.|g, 
m.|.| and m.| 2 . All of the sub-fault identification matrices and M are 
given as follows. 
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E. GENERATION OF X 



The components of X can be obtained by the measurements of the 
normalized |T12l corresponding to the chosen test frequencies. 

Suppose that it is desired to determine which element is faulty in 
a network accessible only at the input and output terminals. The 
diagnosis of the network malfunction can be accomplished using the fault 
identification matrix. 

The value of X will fall in one of the defined fault identification 
regions. 

F. GENERATION OF Z AND FAULT IDENTIFICATION 

The generation of T and the fault identification are best seen by 
the following three cases for example 4.1. 

Case 1. Let the measured values of the components of X at the 
selected test frequencies w.| , ^2 ^nd w^ be 




Z is generated as follows. 

T = X'M 

= [ .6 .7 .9 1. J M 
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Case 2. If the measured values of the components of X at the 
selected test frequencies are 
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Case 3. If the measured values of the components of X at the 
selected test frequencies are 

= 2 . 

^2 = 1.1 

= .98 

Z' = X'M 



= [ 2 . 


1.1 


.98 1 . 


] M 








= [ - 4.34 


4.34 


- 3.7 
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4.54 J 


= t ^1 


^2 
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^5 
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^8 


^9 


^10 


^11 


^12 ] 



From the defined fault identification regions, Table 4-3a, the 
following facts are known for each case. 

Case 1. The algebraically largest component in T is the component z.|. 
The network- is operating in R-1 which means Eal is deviated less than 
-20 percent from the nominal value. In example 4.1, R.| is out of toler- 
ance. 

Case 2. The component Zg in T is the algebraically largest component. 
The netv/ork is operating in R-6. The element Ea2 is deviated more than 
60 percent from the nominal value. In example 4.1, L is out of tolerance. 

Case 3. The component Zg in T is the algebraically largest component 
again but the noticeable fact is that the value of the component Zg is 
different from that of the previous case. The relative value of the 
component of T may suggest the degree of violation of the faulty element. 
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The method discussed in the foregoing illustrations is called the 
single element fault identification method. A combined element fault 
identification, abbreviated by C. E. F. T. , method is a method in v/hich 
more than one element fault is identified. 

The algorithm of the single element fault identification method is 
suitable for computer implementation. If the fault identification 
regions are properly defined, then the construction of M and the genera- 
tion of Z are matrix operations. The fault identification is the 
comparison of the values of the components in T. 

A program for the single element fault identification, designated by 
program 1, is v/ritten as shown in Appendix A and its implementation is 
annexed in Appendix B. The results of the computer applications for 
example 4.1 are. shown in Fig. 4-3, Fig. 4-4 and Fig. 4-5. The same 
procedures are used for example 4.2 and the results of the computer 
applications are shown in Fig. 4-6, Fig. 4-7 and Fig. 4-8. 
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DEFINED BGUNDARY PCINT VALUES 
NUMBER CF FREQUENCIES NV= 3 
NUMBER OF REGIONS NR= 12 



LOWER B.P.V. 

0.0 

0.0 

0.0 

1.20000 

1.20000 

1.06000 

■ 0.0 

0.80000 

1 .02000 

0.0 

0.0 

1.02000 

1.10000 

1.C5000 

0.92000 

1. 20000 

1.05000 
0.92000 

0.0 

0.0 

0.92000 

1.20000 

1.20000 

1.02000 

0 . eoooo 
0.80000 
0.0 

0.0 

0.80000 

0.0 

1.10000 

1.05000 

1.06000 

1.20000 

1.05000 

1 .06000 



UPPER B.P.V. 

0.80000 

0.80000 

0.92000 

100.00000 

100.00000 

100.00000 

0.80000 
0.95000 

1 . 06000 

0.80000 

0.80000 

1.06000 

1.20000 

1.20000 

0.93000 

100.00000 

1.20000 

0.98000 

0.80000 

0.80000 

0.98000 

1 00.00000 

100.00000 

1 . 06000 

■ 0.90000 
0.95000 
0.92000 

0.80000 
0. 95000 
0. 92000 

1.20000 

1.20000 

100.00000 

100. 00000 

1.20000 

100.00000 



Fig. 4-3. 

Fault identification regions 
printed by computer for Example 4.1 
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FAULT IDENTIFICATION MATRIX 
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generated by computer for Example 4.1. 
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DEFINED BOUNDARY POINT VALUES 
NUMBER OF FREQUENCIES NV= 3 
NUMBER OF REGIONS 



LOWER B.P.V. 

0.0 

0.93000 

0.77000 

0.0 

0.82000 

0.77000 

0.0 

0.0 

0.77000 

1.35000 

1.02000 

1.01000 

0.70000 

0.82000 

0.0 

0.70000 

0.0 

0.0 

0.0 

0.0 

0. C 

1. C3000 
1.05000 

1.30000 

1.35000 

0.82000 

0.77000 

1.35000 

0.0 

0.77000 

0.0 

1.20000 

1.01000 

0.95000 

1.20000 

1.30000 

0.95000 

0.0 

0.0 



NR= 13 



UPPER B.P.V. 

0.70000 

0.95000 

0.99000 

0.70000 

0.93000 

0.99000 

0.70000 

0.82000 

0.99000 

100.00000 

1 . 05000 

1.30000 

0.95000 

0.93000 

0.77000 

0.95000 

0.32000 

0.77000 

0.70000 

0.82000 

0.77000 

1.35000 

1.30000 

100.00000 

100.00000 

0.93000 

0.99000 

100.00000 

0.82000 

0.99000' 

0.70000 

100.00000 

1.30000 

0.98000 

100.00000 

100. 00000 

0.98000 

0.82000 

0.77000 



Fig. 4-6. 



Fault identification regions 
printed by computer for Example 4.2. 



FAULT IDENTIFICATION MATRIX 
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V. COMBINED ELEMENT FAULT IDENTIFICATION 



Mul ti -element fault identification, abbreviated by M.E.F.I., is a 
somewhat complicated task. However, if multi-element fault simulation, 
abbreviated by M.E.F.S., for a certain system can be made, M.E.F.I. is 
possible. 

Consider the number of times of M.E.F.S. necessary for a M.E.F.I.. 
For a 3-element system, using the abbreviation for the combination of m 
elements out of n elements, nCm, the number of times to be simulated is 
as follows. 



3C1 = 3 — one element fault combinations 
3C2 = 3 — two elements fault combinations 
3C3 = 1 — three elements fault combination 

Including every element deviation of 10 times, each element deviation 
from -90 percent to +90 percent of the nominal value with 20 percent of 
step size, the number of fault simulations necessary for a M.E.F.I. is 
as many as 



3C1 


10 


= 


30 


for 


one element fault simulations 


3C2 


10 


10 


300 


for 


two elements fault simulations 


3C3 


10 


10 10 = 


1000 


for 


three elements fault simulations 



1330 for total number of fault simulations 

Thus, a total of 1330 mul ti -element fault simulations are required even 
in a simple 3-element system. 

The number of possible M.E.F.S. required for a M.E.F.I. is shown in 
Table 5-1 and it becomes astronomical as the number of elements increases. 
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No. of element 
2 

3 

4 

5 

6 



No. of M.E.F.S. 

no 

1330 

14640 

161050 

1771560 



Table 5-1. 



This astronomical increase of M.E.F.S. discourage^ one to make any 
further investigation. However, another method can take its place as is 
explained as follows. A complex system can be partitioned properly such 
that the probability of failure involving more than one element at a 
time is usually, small . 

In this section, the combined element fault identification (C.E.F.I.) 
method is introduced. A complex network shown in Fig. 5-1 is partitioned 
properly into two subnetworks shown in Fig. 5-2. Two subnetworks NW-1 
and NW-2 are regarded as independent and the single element fault identi- 
fication method can be applied to each netv/ork. 




Fig. 5-1. 



61 




Fig. 5-2. 



Suppose that the fault identification matrices for NW-1 and NW-2 are 
constructed and designated by 



Ml 


for 


NW-1 




M2 


for 


NW-2 




Let M be defined as 










M = 


'mi I 


0 ' 




(5-1) 




1 

0 . 


M2. 




The measured values of the components 


of X can 


be obtained with 


respect to a set of selected test frequencies for 


NW-1 and NW-2 independ- 


ently and designated as follows. 










^ir ^12’ ^13’ ’ ^In’ 




for 


NW-1 


^21’ ^22’ ^23’ ’ ^2m’ 




for 


NW-2 



where n and m are the number of selected test frequencies for NW-1 and 
NW-2 respectively. Then Z are generated as follows. 
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r 



V = X'M = XT , X2‘ [Ml ' 0 



XT Ml, X2'M2 



0 1 M2 



^ ^2 ^3 ^rd ^(rd+ 1 ) ^ 



where 



P = XT, X2 



^ ^11 ^12 ^13 






ZT = XT Ml = [ Z 3 





r is the total number of fault identification regions for both networks 
and rd and r-rd are the number of fault identification regions for NW-1 
and NW-2 respectively. Z will indicate two faulty elements. One from 
Z1 indicates a fault in NVI-1 and the other one from 12 indicates a fault 
in NW-2. 

This is the C.E.F.I. method for two elements faul t .identification. 

The C.E.F.I. method can be extended to isolate multi -element faults. 

For example, suppose that the network shown in Fig. 4-la is replaced 
by NW-1 and the network shown in Fig. 4- lb is replaced by NW-2 for 
C.E.F.I. and the fault identification matrix for NW-1 is Ml and the 
fault identification matrix for NW-2 is M2. The resulting M is constructed 
by Eq. (5-1) and shown in Fig. 5-3. 

Let the measured values of the components of X be 



63 



XI' 


[ Xii x ^2 


Xi 3 1. 


] = 


I 0.6 


0.7 


0.9 


1. 


] 


X2' 


[ X 2 ^ X 22 


X23 1 . 


] = 


[ 0.6 


0.94 0.8 


1. 


] 


X' 


= XI', X2' = 


I 0.6 


0.7 


0.9 


1. 


0.6 


0.94 


0.8 1. ] 



V becomes 

Z' = X'M = XI 'Ml, X2'M2 



[ 3.18 


- 3.18 


2.54 


2.74 


- 0.26 


- 1.46 


3.14 


- 2.86 


2.58 


2.98 


- 0.98 


- 2.18 


3.89 


3.87 


3.63 


- 1.11 


3.61 


3.37 


3.57 


- 0.83 


- 0.15 


- 0.39 


1.45 


- 0.45 


2.67 ] 
















[ Zi • Z2 


^3 


^12 


^13 


“ ^25 


] 







One fault identification region from R-1 R-12, represented by 

2 ^ 2^2 respectively, indicates a fault in NW-1 and the other one 

fault identification region from R-1 3 R-25, represented by 2^^ "" 

— 225 respectively, indicates a fault in NV/-2. In this example, the 

albegraical ly largest component among 2 ^ 2^2 z.| and the largest 

component among 2 -j 2 Z 25 is z^^- The fault identification region 

R-1 indicates that Eal in NW-1 is deviated less than -20 percent from 
the nominal value. The fault identification region R-1 3 represents R-1 
in NW-2 because NW-1 has 12 fault identification regions and indicates 
that Ebl is deviated less than -28 percent from the nominal value. Two 
more examples are presented in Fig. 5-4. 

A program for C.E.F.I., designated by program 2, is written as shown 
in Appendix A and its implementation is annexed in Appendix B. 
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VECTOR X VECTOR 
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Table of vector-Z for C.E.F.I. 



VI. EXPERIMENTS 



A. PASSIVE LINEAR NETWORK EXPERIMENT 

The network shown in Fig. 6-1 (identical to that of Fig. 4-la) v/as 
used in the experiment. 



1 



a 



I'o 



L 




o 2 



■ 02 ' 



R.J = 2 ohms 
L = 2 henry 

C = 3 farads 

R 2 = 4 ohms 



Fig. 6-1. 

Passive network 

The following relationships are used for frequency and magnitude 
seal ing: 



R* = b R 
L* = b/a L 
C* = 1/ab C 
f* = a f 

where "a" is the frequency scale factor (dimensionless, real and 
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positive constant), "b" is the magnitude scale factor (dimensionless, 
real and positive constant) and a scaled nominal element value is 
identified by an asterisk superscript. 

With a = 10000 and b = 100 , the scaled nominal element values 

were 

1/Eal = = 200 ohms 

1/Ea2 = L = 20 mill i henry 
Ea3 = C = 3 micro farads 

1/Ea4 = R 2 “ 400 ohms 

and a set of 3 selected test frequencies became 

w-| = 20000 radians/sec, f.j = 3185 Hz 

W 2 = 10200 radians/sec, f^ = 1624 Hz 

= 200 radians/sec, f^ = 31.8 Hz 

The experimental results for fault simulations at a set of selected 
frequencies were tabulated in Table 6-1 for the network shown in Fig. 
6-1. Comparing Table 6-1 with Table 4-la and taking into consideration 
the accuracy of the measuring devices, it is concluded that the experi- 
mental results in Table 6-1 and the theoretical results in Table 4-la 
are close. 

Consider the following four cases. 

Case 1. If the resistor, was changed to 130 ohms from the 
nominal value of 200 ohms, the magnitudes of T12 were 0.0935 , 0.222 

and 0.755 at a set of 3 selected test frequencies f.j , f 2 and f^ 
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respectively. Since the nominal magnitudes of T12 were the same as used 
in Section IV-C, the normalized magnitudes of T12 became 

= 1.562 

X 2 = 1.515 

X 3 = 1.138 

Case 2. If the inductor, L, was changed to 30 milli henry from the 
nominal value of 20 milli henry, the magnitudes of T12 at a set of 3 
selected frequencies were 0.0462 , 0.131 and 0.664. The normalized 
magnitudes of T12 were 

■x^ = 0.777 

X 2 = 0.891 

X3 = 1.0 

Case 3. If the capacitor, C, was changed to 4 micro farads, the 
magnitudes of T12 were 0.448 , 0.109 and 0.664 and the normalized 
magnitudes of T12 were 

x^ = 0.748 
X 2 = 0.743 
X 3 = 0.98 

Case 4. If it was made that the resistor, 839 changed to 250 ohms, 
the magnitudes of T12 were 0.0453 , 0.13 and 0.56 and the normalized 

magnitudes of T12 were 
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0.757 



X2 = 0.889 
X3 = 0.84 

Using tfie fault identification matrix. Fig. 4-4, Z and the fault 
identification were generated as shwon in Fig. 6-1. 

From Fig. 6-1 and the defined fault identification regions. Table 4-3a, 
the following facts were known for each case. 

Case 1. The component was the largest component in Z and X 
belonged to the fault identification region R-2. Eal was deviated more 
than 20 percent from the nominal value. R.j was out of tolerance. 

Case 2. The network was operating in R-3. Ea2 was deviated to 
somewhere between -30 percent and -50 percent from the nominal value. 

L was out of tolerance. 

Case 3. The network was operating in R-7 and Ea3 was deviated more 
than 25 percent. C was out of tolerance. 

Case 4. Ea4 was deviated more than 50 percent which was indicated 
by R-10. 

B. ACTIVE LINEAR NETWORK EXPERIMENT 

The network shown in Fig. 6-2 was used for active linear network 
experiment. The frequency response of the active network was drawn in 
Fig. 6-3 using the computer program ECAP. A set of selected test fre- 
quencies was chosen as follows. 

f^ = 10 Hz 

f2 = 100 Hz 

f3 = 1000 Hz 
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Experimental fault simulated data for Passive network 
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Table of Vector-Z for passive network experiment. 
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(b) small signal ECAP 

equivalent circuit for Q1 



= 100 K ohms 

R 2 = 10 K ohms 

R^ = IK ohms 

R^ = 5 K ohms 

= 5 pf 

C 2 = 50 pf 

Q1 — 2N1414 

H.g = 1300 ohms 

1/H^g = 36 K ohms 
= 2.9 10"^ 

Hfe= 42 

Coe =26 Pf 
614 - 



Fig. 6-2. 
Active network 



Since there was no computer program available for fault simulations 
of the active case, the necessary data were obtained from the experi- 
mental results. The fault simulated data were tabulated in Table 6-2. 
The fault simulation curves were drawn in Graph 6-1, Graph 6-2 and Graph 
6-3. The normal regions were defined in Table 6-3 and the fault identi- 
fication regions were defined in Table 6-4. The fault identification 
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matrix was generated by the program 1 (in Appendix A) and is shown in 
Fig. 6-5. 

Z and the fault identification were obtained for the following 
three cases. 

Case 1. Let the resistor, R 2 , be changed to 14 K ohms from the 
nominal value of 10 K ohms. The measured magnitudes of T12 were 11.71 , 

91.0 and 179.0 at 3 test frequencies f-j , f 2 and f^ respectively. 

Then the normalized magnitudes of T12 became 

X.J = 1.03 

X 2 = 1.08 

X 3 = 1.26 

Case 2. If the resistor, R^, was changed to 7 K ohms from the 
nominal value of 5 K ohms, the measured magnitudes of T12 were 15.9 , 

126.0 and 190.2 at 3 corresponding test frequencies. The normalized 
magnitudes of T12 were 

X., = 1.39 

X 2 = 1.5 
X 3 = 1.34 

Case 3. If the capacitor, C 2 , was changed to 35 micro farads from 
the nominal value of 50 micro farads, the measured magnitudes of T12 
were 9.5 , 75.0 and 143.0 at 3 corresponding test frequencies. The 
normalized magnitudes of T12 were 
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0.834 



0.892 
1.007 ' 

Using M, shown in Fig. 6-5, Z and the fault identification generated 
by the computer (program 1) are shwon in Fig. 6-6. 

From the defined fault identification regions, Table 6-4, the follow- 
ing facts were known for each case. 

Case 1. The component was the algebraically largest component in 
Z. R-3 indicated that R 2 was deviated more than 30 percent. 

Case 2. The fault identification region R-8 indicated that R^ was 
deviated more than 20 percent. 

Case 3. The fault identification region R-13 indicated that C 2 was 
deviated between 40 percent and 70 percent. 

C. PASSIVE AND ACTIVE NETWORK C.E.F.I. EXPERIMENT 

Suppose that the networks shown in Fig. 6-1 and Fig. 6-2 v/ere 
considered as an interconnection of two-port network consisting of two 
cascade connected subnetworks and regarding Fig. 6-1 as NW-1 and Fig. 6-2 
as NW-2, The fault identification matrix for C.E.F.I. experiment were 
shown in Fig. 6-7. 

Applying the same procedures as in the C.E.F.I. method previously 
discussed, the netv/ork malfunctions could be isolated. The C.E.F.I. was 
accomplished for three cases with the same X used in the passive and 
active netv/ork single element fault identification experiments. The 
resulted fault identification was shown in Fig. 6-8. 
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NIV0 
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Frequency response for active network. 



Element DEV Element MAG-Tl MAG-T2 MAG-T3 NCR-Tl N0R-T2 N0R-T3 

(percent) Value (at f^) (at f^) (at f^) 
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Experimental fault simulated data for Active network. 
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Experimental fault simulation curve for Active network. 
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Experimental fault simulation curve for Active network. 
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Graph 6-3. 

Experimental fault simulation curve for Active network. 
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FAULT IDENTIFICATION REGIONS 
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DEFINED BOUNDARY POINT VALUES 
NUMBER OF FREQUENCIES NV= 3 
NUMBER OF REGIONS NR= 14 



LOWER B.P.V. 

0.80000 

0.0 

0.0 

0.0 

0.0 

0.0 

1.02000 

1.07000 

1.20000 

1.05000 

1.07000 

1.20000 

1.20000 

1.07000 

1.20000 

0.92000 
0. 80000 
0.66000 

0.0 

0.0 

0.66000 

1.20000 

1.20000 

1.20000 

1.05000 

1.07000 

0.84000 

1.20000 

1.20000 

0.84000 

1.20000 

1.20000 

0.66000 

0.92000 

0.93000 

1.01000 

0.80000 

0.80000 

1.01000 

0.0 

0.0 

1.01000 



UPPER B.P.V. 

0.92000 

0.80000 

0.66000 

0.80000 

0.80000 

0.66000 

1.05000 

1.20000 

100.00000 

1.20000 

1.20000 

100.00000 

100.00000 

1.20000 

100.00000 

0.97000 
0.93000 
0. 84000 

0.80000 
0. 80000 
0.84000 

100.00000 

100.00000 

100.00000 

1.20000 

1.20000 

0.99000 

100.00000 

100.00000 

0.99000' 

100.00000 

100.00000 

0.84000 

0.97000 

0.97000 

1.20000 

0.92000 

0.93000 

1.20000 

0. 80000 
0.80000 
1.20000 



Fig. 6-4. 

Fault identification regions 
printed by computer for Active netv/ork. 
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FAULT IDENTIFICATION MATRIX 
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VECTOR X VECTOR 
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Table of vector-Z for Active network. 
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Table of vector-Z for C.E.F.I. experiment. 



i 




VII. FEASIBILITY OF FAULT PREDICTION 



The feasibility study of fault prediction is quite complicated. 
However a brief discussion of this topic is included in this section. 
Suppose that M of a certain system is 
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... 



and X' is 

X “ E ^1 ^2 ^3 ~ ^ 

7' becomes 



V = X'M = [ X.J X 2 x^ 1 . ] M 



= t 2l 22 — ^ ] 



where 



^0 



^l^lj 



^2^2j 



n nj 



A. 

J 



j=l ,2,— ,r 



X. is positive real number ( i=l,2, — ,n) and a., and A. are real 
1 1 j j 

constants. 
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Since z. is linear functions of x.s, the variation of the value of 

J * 

z. is monotonic. 

J 

If the system is faulty and z. is algebraically the largest component, 

J 

the system is operating in R-j and an element which is out of tolerance 
can be isolated from the predefined fault identification regions. 

Taking the network shown in Fig. 4-la as an example for the fault 
prediction, there are 8 boundary point values between the normal region 
and the faulty region at a test frequency. The largest components in Z 
at these boundary point values are shown in Table 7-1. From this table 

Element x, ^2 ^3 TheJIargest component 

(percentage) boundary) ^ boundary 



Eal 


(-20) 


(0.80 


0.80 


0.91) 


1.90 


Eal 


(+20) 


(1.20 


1.20 


1.06) 


1.86 


Ea2 


(-30) 


(0.80 


0.95 


1.02) 


1.48 


Ea2 


(+25) 


(1.10 


1.05 


0.98) 


1.28 


Ea3 


(-18) 


(1.20 


1.20 


1.02) 


1 .78 


Ea3 


(+25) 


(0.80 


0.80 


0.98) 


1.78 


Ea4 


(-20) 


(1.10 


1.05 


1.06) 


1.36 


Ea4 


(+25) 


(0.90 


0.95 


0.92) 


1 .40 



Table 7-1. 

the following observations are made: 

1) If the deviation of Eal element is to be decreased less than 
-20 percent, the largest component in T must be algebraically greater 
than 1.90-, if the deviation of Eal element increased more than +20 
percent, the largest component in Z must be greater than 1.86. The same 
idea can be applied to all elements. 
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2) If z. for all j is smaller than 1.28, the smallest value in 
Table 7-1, the network is working in normal condition. 

3) If z. for all j is larger than 1.90, the largest value in Table 

J 

7-1, the network is working in one of the faulty regions. 

4) If the largest component z. in T is in between 1.28 and 1.90, 

• J 

the network may or may not be in a normal region. For instance, if the 

largest component z. is 1.50 and the network is still in normal condition, 

J 

then it is known that Ea2 and Ea4 are not faulty and Eal or Ea3 may begin 

to approach one of the faulty regions. If the largest component z. is 

J 

1.50 and the network is in faulty condition, then it is known that Ea2 
is faulty. 
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VIII. CONCLUSION 



A fault identification technique has been investigated in this paper. 
Although the examples used are simple, the fault identification proce- 
dures presented here may be extended to more complex systems. 

The theoretical basis of the proposed method is described and the 
technique has been applied to both passive and active networks. 

The main differences in the approach used in this. paper for fault 
identification over the gain signature approach cited in References [3], 
[8] and [10] are summarized as follows. 

1) . Fault identification regions are defined uniquely according to 
the fault simulation curves which are constructed from the system sensi- 
tivity instead of assignment of gain signature. 

2) . Direct application of the measured values for fault isolation 
can be accomplished instead of gain signature conversion. 

3) . Vector-Z indicates an exact faulty element automatically rather 
than a gain signature comparison. 

4) . The selection of a set of test frequencies is not critical as 
in the gain signature approach. 

5) . By applying the matrix operation to the fault identification 
matrix of the single element fault identification method, one can detect 
not only a fault of the one-element type but also a multi -element fault 
in the given system. 

The proposed method uses only one network function. However, it 
gives no ambiguities as demonstrated in the examples. 
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DIMENSION BI30»50) tC( 50»50) tY(50t30) ,JJ(50) 
READ(5,111) NI,NRtND 

111 F0RMAT(3I10) 

RE AO (5, 112) ( ( B( II » J) , I 1 = 1 ,NI ) , J = 1,NR) 

112 FORMAT ( 8F10.0) 
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DO 172 11=1, NI 

17 2 C(KK,J)=C(KK,J)+Y(KK,II)’!'B{ II,J 
IF(KK-1 ) 173,173,174 
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581 DO 582 M=1,ND 

582 JJ(M)=0 
U=C( KK, J) 
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APPENDIX B 



IMPLEMENTATION OF THE COMPUTER PROGRAMS 



INPUT DATA 
PROGRAM 1 



Card 

number 


Card 

column 


Information 


Form 


1 


1 


— 10 


NV number of selected frequencies 


i nteger 




11 


— 20 


NR number of fault identification 
regions 


II 


2 


1 


— 10 


lower B.P.V. of x-| for fault 
identification region R-1 


decimal 




11 


- 20 


upper B.P.V. of x-j for fault 
identification region R-1 


II 


3 


1 


- 10 


lower " X 2 " 


II 




11 


- 20 


upper " X 2 " 


II 


NV+1 


1 


— 10 


lower " x^y " 


II 




11 


— 20 


uPPe^' " Xnv " 


II 


NV+2 


1 


- 10 


lower " x-| " 

identification region R-2 


II 




11 


- 20 


upper " X] " 

identification region R-2 


II 


2NV+1 


1 


- 10 


lower " x^y " 


II 




11 


— 20 


upper " x^^ " 


II 
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2NV+2 


1 


— 10 


lower " x*| 

identification 


II 

region 


R-3 


II 




11 


- 20 


upper " x-j 

identification 


II 

region 


R-3 


II 



NR-NV+1 

NR-NV+2 


1 


_ _ 


10 


NM number of measurements 


integer 


NR-NV+3 


1 


— 


10 


magnitude value of x-j 


decimal 




11 


— 


20 


II Y 

^2 


II 




21 


— 


30 


II Y 

^3 


II 


NR-NV+4 


1 


NV 


10 


1. 

magnitude value of x^ 


II 




11 


— 


20 


" ■ ^2 


II 






NV 




1. 





NR-NV+2+NM 



99 



PROGRAM 2 



Card Card 

number column 



Information 



Form 



1 


1 


-- 


10 


NI 


number of M 


rows 


integer 




11 


— 


20 


NR 


number of M 


columns 


II 


2 


1 


— 


10 


entry of col-1 , 


row-1 of M 


decimal 




11 


— 


20 




II 


row-2 " 


II 




21 


— 


30 




II 


row- 3 " 


II 














row-ni " 


II 


3 


1 


— 


10 


entry of col -2, 


row-1 of M 


II 




11 


— 


20 




H 


row- 2 " 


II 














row-ni 


II 


NR+1 










• 






NR+2 


1 


— 


10 


NM 


number of measurements 


integer 


NR+3 


1 


— 


10 


magnitude value 


of X-| 


decimal 




11 


— 


20 




II 


^2 


II 






NV 






1. 






NR+4 


1 


— 


10 


magnitude value 


of X-| 


II 




11 


— 


20 




II 


^2 


II 



NV 1. 



NR+2+NM 

* Open end value may be used +100. for upper boundary and 0, for 
lower boundary. 

* Appropriate modification may be done. 
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